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ABSTRACT 

The  feasibility  of  pulse  torquing  a  pendulous  accelerometer 
in  a  ternary  mode  is  investigated.     A  comparison  is  made  between 
basic  ternary  and  binary  systems.     It  is  shown  that  a  ternary  sys- 
tem is  basically  more  stable  than  a  binary  system  in  that  stable 
oscillations  are  smaller.     Thus  errors  in  velocity  indication  are 
less  for  the  ternary  system. 

The  describing  equations  of  a  pulse  torqued  pendulous  ac- 
celerometer are  transformed  into  a  non-dimensional  form.     By  so 
doing,   general  conclusions  are  reached.     Phase  plane  analysis  is 
used  to  investigate  the  dynamic  performance  of  the  device.      An 
error  analysis  based  upon  the  phase  plane  analysis  follows. 

After  conditions  necessary  for  larger  limit  cycles  are  inves- 
tigated, a  method  is  devised  for  selecting  system  parameters  such 
that  larger  limit  cycles  are  prevented. 

Finally,  methods  of  compensation  and  error  correction  using 
switching  logic  are  investigated.     Representative  switching  logics 
are  presented. 

Thesis  Supervisor:    W.    E.   Vander  Velde 

Title:  Assistant  Professor  of 
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CHAPTER  1 
INTRODUCTION 

1. 1.     HISTORY 

A  basic  requirement  of  any  inertial  navigation  system  is  that  the 
system  be  able  to  obtain  velocity.     One  common  device  developed  to 
obtain  velocity  inertially,   is  the  pendulous  accelerometer.     In  such  a 
device,   a  torque  generator  current  required  to  keep  a  pendulum  at 
null,   is  proportional  to  acceleration  along  a  sensitive  axis.     When  such 
current  is  integrated  with  respect  to  time,   an  analog  of  velocity   is  ob- 
tained. 

With  the  addition  of  a  digital  computer  to  the  inertial  navigation 
system,   a  requirement  is  generated  for  an  analog  to  digital  converter. 
Such  converters  are  inherently  limited  in  accuracy  of  conversion. 

As  an  alternative  to  providing  a  converter,   a  binary  pulse  torqued 
integrating  accelerometer  has  been    under  study  and  development.    Such 
an  accelerometer  has  an  output  which  consists  of  a  train  of  pulses. 
These  pulses  may  be  fed  directly  to  a  digital  computer.     Each  pulse 
represents  an  incremental  velocity  change  along  the  sensitive  axis  of  the 
accelerometer.     By  summing  the  train  of  pulses,    taking  into  considera- 
tion their  algebraic  sign,   a  velocity  indication  is  obtained. 

In  the  operation  of  a  binary  pulse  torqued  integrating  accelerom- 
eter,   certain  errors  are  introduced  which  are  caused  by  inherent  lag 
in  the  pendulum  servo  system. 

1.  2.     A  PENDULOUS  ACCELEROMETER 

In  general,   a  pendulous  accelerometer  consists  of  a  pendulous 
mass  suspended  in  a  viscous  fluid  (Figure  1-1).     The  pendulous  mass  is 
constrained  to  rotate  about  a  single  axis.     An  error  detector  is  used  to 
measure  the  deviation  of  the  pendulous  mass  from  a  reference  axis. 
After  amplification,   the  error  signal  is  used  to  command  a  torque  gen- 
erator in  such  a  way  that  the  pendulum  is  driven  back  toward  its  null 
position.     The  applied  torque  may,   or  may  not  be,   a  linear  function  of 
null  angle. 
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The  equation  of  motion  of  the  pendulum  is 

2 
HA  HA 

J  iJL  +  c  —   =  M.     -  a  P  cos  A  +  a„  .    P  sin  A  (1.1) 

2  tg  RA 

at 

where 

P  =  Pendulosity  of   the  pendulum 

A  =  Angular  position  of  the  pendulum 

J  =  Moment  of  inertia  of  the  pendulum 

C  =  Coefficient  of  viscous  friction 

M,      =  Torque  Generator  torque 
tg 

a  =  Component  of  acceleration  along  the  input  axis  of  the 
accelerometer 
aR  a    =  Component  of  acceleration  along  the  Reference  Axis 
of  the  accelerometer 
t  =  Time. 
For  the  purpose  of  this  report,   acceleration  is  defined  as  the 
second  derivative  of  position  with  respect  to  time  as  viewed  in  a  free- 
fall  non-rotating  coordinate  system.     This  is  not  a  restriction  in  that 
Equation  (1.  2)  is  valid  in  any  coordinate  system  if  the  acceleration,   a, 
is  defined  to  include  gravitational  attraction  and  apparent  accelerations 
associated  with  rotating  frames. 

The  null  angle,   A,   is  maintained  sufficiently  small  to  justify  small 
angle  approximations  being  made.     With  this  in  mind,   Equation  (1. 1) 
becomes: 

2 

j5L4  +  C—   =Mtg-aP  +  aRAPA  <L2> 

dr  dt 

Rearranging  and  integrating: 

v=\      adt=-      f     M      dt-—  —  -—  A  +  \      a         Adt  (1.3) 

Jo  P    Jo       tg  P  dt         P  Jo      K 

where  v  is  the  velocity  increase     from  time  zero  to  time  t.     dA/dt  and 
A  are  assumed  zero  at  time  equal  to  zero. 
By  definition: 
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1  ft 

v.  =  —    \       MA     dt  =  indicated  velocity  (1.  4) 

i      p  J0         tg  J 


e     =  —  — —  +  —  A  =  null  error  (1.  5) 

n 

P  dt         P 


■-r 


e     =  -   \       aTDA  A  dt  =  cross-coupling  error  (1.6) 


o 


e  =  e     +  e     =  total  error.  (1.7) 

xi         c 

Substituting  (1.5)  and  (1.6)  into  (1.3): 

v=—    \M,     dt-e-e.  (1.8) 

p    J0       tg  n         c 

1.  3.     A.  BINARY  PULSED  INTEGRATING  ACCELEROMETER 

Figure  1-2  is  a  block  diagram  of  a  binary  torqued  integrating  ac- 
celerometer.     The  pendulum  is  always  torqued  by  a  current  of  constant 
magnitude.     The  sense  or  direction  of  the  current  through  the  torque 
generator  may  be  changed,   however.     The  error  detector  is  sampled 
at  a  rate  synchronous  with  a  computer  clock.     At  each  sample     instant, 
the  current  to  the  torque  generator  is  switched  or  not  switched  as 
necessary  to  drive  the  pendulum     toward  null.     Lag  is  introduced  into 
the  restoring  system  by  the  use  of  a  discrete  time  between  samples. 

The  state  of  the  applied  torque  is  based  upon  the  error  signal  at 
a  sample  instant.     Practical  error  detectors  have  a  finite  insensitive 
zone  near  their  null  position.     If  at  a  sample  instant,   the  pendulum  is  in 
the  insensitive  zone  (hereafter  called  the  dead  zone),   the  applied  torque 
is  maintained  in  a  sense  based  upon  the  last  available  information.    Thus 
additional  lag  is  introduced. 

As  a  result  of  lag  in  the  error  channel,  the  pendulum  oscillates  in 
stable  modes  or  limit  cycles.  These  limit  cycles  if  large  allow  the  null 
error,   e   ,   to  be  several  velocity  increments  in  magnitude. 

XI 
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Figure  1-2.     A  Binary  Pulsed  Pendulous 
Integrating  Accelerometer. 
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The  chapters  to  follow  investigate  the  feasibility  of  reducing  the 
errors  in  a  pulse  torqued  integrating  accelerometer  by  torquing  the  pen- 
dulum in  a  ternary  mode. 
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CHAPTER  2 
A  TERNARY  PULSED  INTEGRATING  ACCELEROMETER 

2.1.  GENERAL 

The  size  of  the  limit  cycle  of  a  pulsed  accelerometer  may  be  re- 
duced by  one  of  two  methods.     First,  the  lag  which  causes  the  ex- 
cessive limit  cycles  may  be  reduced.     Second,  lead  may  be  introduced 
to  compensate  for  lag  already  in  the  system. 

If,   in  a  pulse  torqued  accelerometer,   the  pendulum  null  angle  is 
such  that  the  acceleration  reaction  torque  is   in  a  direction  to  reduce  the 
null  angle,  it  can  be  anticipated  that  the  null  angle  will  eventually  be 
driven  to  zero  without  the  use. of  additional  torque.     Thus,   the  use  of 
applied  torque  to  aid  acceleration  reaction,  coupled  with  the  discrete 
time  increments  at  which  it  must  be  applied,   introduces  unnecessary  lag 
to  the  system.     This  elementary  observation  prompted  this  investigation. 

If  torquing  to  aid  acceleration  reaction  torque  is  to  be  prevented, 
a  "no  torque"  mode  of  operation  must  be  introduced.     Since  both  posi- 
tive and  negative  accelerations  are  to  be  expected,   positive  and  negative 
torque  modes  must  be  retained.     Thus  a  ternary  mode  of  operation  is 
suggested. 

2.2.  COMPUTER  CONSIDERATIONS 

The  computer  most  usually  used  in  inertial  guidance  systems  is 
the  Digital  Differential  Analyzer  (DDA).     The  Digital  Differential  Analyzer 
can  be  mechanized  to  handle  binary  or  ternary  inputs.     The  ternary 
scheme  is  more  accurate   in  that  it  has  less  round-off  errors,  but  re- 
quires more  equipment  to  realize. 

The  DDA  must  accept  information  at  a  clock  rate.     Thus  the  re- 
quirement that  the  torque  pulses  be  synchronous  with  a  computer  clock 
will  be  maintained  in  a  ternary  system. 

2.  3.     INFORMATION  AVAILABLE 

The  accuracy  of  a  pendulous  accelerometer  depends  upon  the  devia- 
tion of  the  pendulum  from  its  null  position  being  maintained  small.     With 
such  small  deviations  allowed,  pendulum  velocity  is  extremely  ha^d   to 
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meo  sure   through  system  noise.     For  this  reason,   it  is  considered  that 
pendulum  velocity  is  not  directly  available. 

Pendulum  position  is  available  if  the  pendulum  is  not  in  the  error 
detector  dead  zone.     If  the  pendulum  is  in  the  dead  zone,    the  fact  that 
it  is  in  the  dead  zone  is  known.     This  fact  may  be  used  in  switching  logic. 

The  past  history  of  applied  torque  is  available  and  may  be  useful 
in  that  it  represents  the  past  motions  of  the  pendulum.     Since  input  ac- 
celeration is  a  slow  changing  function  with  respect  to  pendulum  motion, 
the  past  history  of  pendulum  motion  might  be  useful  in  predicting  future 
motion.     Thus,  lead  might  be  introduced. 

2.  4.     POWER  SUPPLY  CONSIDERATIONS 

If  torque  pulses  are  to  accurately  indicate  the  same  velocity  quan- 
tum with  repeatibility,   the  torque  generator  current  must  be  maintained 
extremely  constant.   'The  current  taken  from  the  power  supply  is  con- 
stant in  a  binary  system.     A  ternary  system  requires  that  the  current 
be  switched  on  and  off.     Because  of  the  variation  in  current  taken  from 
the  power  supply,  a  better  power  supply  would  be  required  for  the  ternary 
system. 

An  alternative  to  switching  the  power  supply  off  might  be  to  switch 
it  to  a  dummy  load.     At  any  rate  a  ternary  system  would  require  more 
or  better  power  supply  equipment. 

2.5.     A  BASIC  TERNARY  PULSED  INTEGRATING  ACCELEROMETER 

Figure  2-1  shows  a  basic  ternary  pulsed  integrating  accelerometer. 
Switching  logic  is  based  upon  pendulum  position.     If  at  a  sample  instant 
(clock  pulse),   the  pendulum  is  not  in  the  error  detector  dead  zone,   the 
pendulum  is  torqued  for  the  next  sample  interval,   in  a  direction  toward 
null.     If  the  pendulum  is  in  the  dead  zone  at  a  sample  instant,   the  pen- 
dulum is  not  torqued  for  the  next  sample  period. 

For  a  pulsed  system,   Equation  (1.  8)   can  be  written: 

n 


where 


v=_    2    5iMts-en-ec-ed  (2.1) 

*    i=0 

t     =  time  between  samples 
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Figure  2-1.     A  Ternary  Pulsed  Pendulous  Integrating  Accelerometer, 


+1  when  torquing  in.  a  positive  direction 
6.  ::  <  0    when  not  torquing 

-1   when  torquing  in  a  negative  direction 


t 


and  e ,  =  —    )    6.Mt     -  —    \ 


p  to  *  p  J° 

The  digital  error  is  described  by  Figure  2-2.     It  is  caused  by  an 
instantaneous  addition  of  velocity  quantum  in  the  computer.     The  actual 
torque- time  pulse  has  taken  place  over  a  finite  time  interval. 
The  total  error  is  now  defined  by: 

e  =  e     +  e     +  e  ,  (2.2) 

n         c         d 

2.6.     NON-DIMENSIONAL  FORM 

It  is  desirable  to  express  the  equations   of  motion  as  free  of  as 
many  system  constants  as  possible.     By  so  doing,   general  conclusions 
can  be  easier  reached.     By  proper  substitutions  of  relations  in  Table  2-1 
into  Equations  (1.2),   (2.1),   and  (2.2): 

i'+i  =  (6  -  r  {  rRA  0)  (2.3) 


n 


-)      6.  -  E     - 
L       i         n 


V  =   ;      6.  -  E     -  E     -  E  ,  (2.4) 

i=0 


E  =  E     +  E     +  E  ,  (2.  5) 

n         c  d  \    •    / 

where  (/>  _  — x~  . 

dT 

Two  system  parameters  must  be  specified  before  an  analysis  of 
system  performance  can  be  made.     Equation  (2.  3)  will  have  to  be  piece- 
wise  integrated  between  samples.     Therefore  T   ,   the  ratio  of  sample 
time  to  pendulum  time  constant,   must  be  specified.     The  non-dimensional 
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Figure  2-2.     Example  of  Digital  Error. 
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TABLE  2-1 
DIMENSIONAL  TO  NON-DIMENSIONAL  TRANSFORMATIONS 

Parameter  Dimensional         Non-Dimensional  Relation 

Notation  Notation 

Angle  A  £ 


Time 


Acceleration  a  r 

Velocity  v  V 

Velocity  Error  e  E 


$-- 

=  AC 

X. 

MJ 

T  = 

:ct 

J 

r  = 

.  aP 

M 

V 

PC 

=  V 

MJ 

E 

PC 
=  e 

MJ 
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deaci  zone,  <£.,  must  be  specified  in  order  to  tell  where  switching  will 
take  place. 

2.7.     PHASE  PLANE  ANALYSIS 

In  general,   the  acceleration  is  essentially  constant  when  viewed 
over  several  sample  intervals.     Thus,   acceleration  is  considered  con- 
stant in  this  analysis.     For  the  purpose  of  analysis,  r  is  assumed  much 
greater  than  rR  »  </>  in  prescribing  the  motion  of  the  pendulum.     This  does 
not  imply  that  the  cross -coupling  error  is  negligible  but  rather  that  the 
dynamic  behavior  of  the  pendulum  is  not  greatly  affected  by  the  cross- 
coupling  acceleration  torque.     This  type  analysis  cannot  be  made  for 
zero  acceleration  along  the  sensitive  axis  in  that  r_  .   is  the  only  forcing 
function  when  the  pendulum  is  not  in  a  torque  mode. 

With  the  r    »    rR.  </>  approximation,  Equation  (2.3)  becomes: 

i"  +  i  =  (6  -  r)  (2.6) 

Since  6  remains  constant  between  torque  switches,   Equation  (2.6) 
may  be  integrated  as  a  first  order  linear  equation  in  <£,  with  results  being 
valid  between  torque  switches: 

i  =  i0  e"T  +(5  -  r)(l  -  e"T)  (2.7) 

where  (/>     is  the  pendulum  velocity  at  t  =  0. 

4>  =  <j>o  +  iQ  (1  -  e"T)  +  (6  -  r)(  t  +  e"T  -  1)  (2.  8) 

where  <j)     is  the  null  angle  at  t    =0. 

It  is  believed  that  a  phase  plane  (</>  versus  $  plot)  study  lends  the 
best  understanding  of  second  order  nonlinear  equations.     It  is  assumed 
that  the  reader  has  a  basic  understanding  of  phase  plane  analysis.    Since 
an  elastic  restraint  term  is  absent  from  Equation  (2.6),   a  phase  plane 
template  can  be  constructed  for  a  given  value  of  (6  -  r).     Such  a  template 
can  be  used  over  and  over  again  by  sliding  it  along  the  </>  axis.     This  is 
possible  because  the  phase  trajectory  is  independent  of  the  relative 
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position  of  the  </>  =  0  axis.     Time  must  be  plotted  as  a  parameter  along 
the  trajectory  in  order  to  tell  when  switching  will  take  place. 

Since  the  null  error  has  been  defined  as  the  sum  of  angle  and  angle 
rate,  it  can  be  found  for  any  point  on  the  phase  plane  by  projecting  the 
point  in  question  onto  an  axis  which  is  perpendicular  to  all 
</>+#=  constant  lines. 

Adding  Equations  (2.7)  and  (2.8): 

E     =  </)     +  6     +  (6  -  r)  t.  (2.9) 

n      ro      ^o 

Equation  (2.  9)  states  that  the  null  error  is  a  linear  function  of  time 
between  torque  switches. 

2.  8.     PREVENTING  APPLIED  TORQUE  FROM  AIDING 
ACCELERATION  REACTION  TORQUE 

It  is  desirable  that  torque  not  be  applied  in  a  direction  that  will  aid 
the  acceleration  reaction  torque.     Appendix  B  shows  that  if  the  non- 
dimensional  dead  zone,  (/>.,  is  greater  than  the  non-dimensional  sample 
time,   T   ,   a  limit  cycle  cannot  exist  for  zero  input  acceleration.     In  gen- 
eral, if  </>.    is  greater  than  T  ,   the  pendulum  in  steady  state,  will  re- 
main out  of  the  switch  zone  that  will  switch  torque  to  aid  acceleration 
reaction  (Figures  2-3  through  2-12).     For  accelerations  greater  than 
zero,   the  forcing  function  of  Equation  (2.6)   tends  to  bend  the  phase  tra- 
jectory away  from  the  undesirable  zone. 

2.9.     EXAMPLE 

In  any  analysis  of  a  pulsed  accelerometer,   the  non-dimensional 
sample  time  and  the  non-dimensional  dead  zone  are  functions  of  system 
parameters.     For  an  illustrative  example  of  analysis  by  phase  plane 
methods,   the  non-dimensional  sample  time  and  the  non-dimensional  dead 
zone  are  both  taken  as  unity.     The  non-dimensional  dead  zone  is  assumed 
greater  than  the  non-dimensional  sample  time  by  an  infinitesimal  amount, 
such  that  limit  cycles  cannot  exist  for  zero  acceleration. 

Figures  2-3  through  2-12  are  phase  plane  plots  of  the  system  re- 
sponse for  various  values  of  non-dimensional  acceleration.     For  magni- 
tudes of  non-dimensional  acceleration  equal  to  or  greater  than  one,   the 
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restoring  torque  cannot  overcome  the  acceleration  reaction  torque.     The 
system  becomes  unconditionally  unstable  under  such  conditions.    For  ac- 
celerations within  the  designed  range  of  the  accelerometer,   the  non- 
dimensional  acceleration  will  never  have  a  magnitude  greater  than  one. 
Phase  plane  plots  of  a  binary  torqued  system  are  shown  in  Figure 
2-3  through  2-12  for  comparison  purposes.     The  binary  and  ternary  sys- 
tems have  identical  system  parameters.     As  could  be  expected,  the  limit 
cycles  of  the  ternary  system  are  smaller  than  those  of  the  binary  system 
for  the  same  input  accelerations.     Higher  order  limit  cycles  exist  for 
both  systems.     Such  limit  cycles  give  rise  to  errors  in  velocity  indication. 
For  the  associated  accelerations,   all  possible  limit  cycles  have  been  shown 
in  Figures  2-3  through  2-12.     The  conditions  under  which  higher  order 
limit  cycles  may  exist  in  a  ternary  system  are  investigated  in  Chapter  3 
and  Appendix  C. 

2.10.     ERROR  ANALYSIS 

From  phase  plane  studies,  it  is  noted  that  limit  cycles  for  a  given 
input  acceleration  have  a  certain  freedom  in  position  about  the  switching 
line.     The  exact  amount  of  freedom  that  a  simple  limit  cycle  possesses 
has  been  determined  and  plotted  as  a  function  of  T    in  Appendix  C. 

Figures  2-13  through  2-22  are  plots  of  null  error  versus  time  for 
the  phase  plane  plots  of  Figures  2-3  through  2-12.     Again,  a  plot  is  shown 
for  a  binary  system  for  comparison.     The  error  plots  are  not  unique  in 
that  the  initial  conditions  determine  the  exact  location  of  the  limit  cycle. 
The  null  error  is  bounded  by  a  region  determined  by  the  freedom  in  po- 
sition of  limit  cycle,  </>f,  plotted  in  Appendix  C. 

it  is  noted  that  the  peak  null  error  is  always  less  for  the  ternary 
system,  usually  by  as  much  as  one  to  two  velocity  counts. 

The  cross-coupling  error  is  equal  to  the  product  of  the  component 
of  acceleration  along  the  reference  axis,  and  the  average  value  of  the  null 
angle.     Both  systems  are  biased  in  that  the  average  null  angle  is  not  zero 
for  accelerations  greater  than  zero.     In  general,   the  average  null  angle 
is  greater  for  the  ternary  system  if  the  acceleration  along  the  input  axis 
is  small  and  is  greater  for  the  binary  system  if  the  acceleration  along 
the  input  axis  is  large.     The  bias  in  null  angle  is  equivalent  to   the  forced 
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dynamic  error  in  a  pendulous  accelerometer  with  linear  feedback. 

The  maximum  magnitude  of  the  average  null  angle  is  always  less 
than  the  magnitude  of  the  dead  zone.  The  maximum  value  of  the  cross 
coupling  error  is 


c 


SrRA*A  <2-l°> 

max 


In  general,   the  digital  error,    E ,,  is  of  opposite  sign  than  the 
null  error,   thereby  reducing  the  total  error.     The  digital  error  is 
the  error  introduced  by  an  incremental  summation  being  made  in  the 
counter,   instead  of  a  true  integration  of  torque. 

The  error  equations   used  in  this  analysis  have  assumed  that  the 
pendulum  is  at  rest  in  its  null  position  when  the  problem  commences. 
An  exact  analysis  will  require  that  the  error  associated  with  initial 
conditions  be  subtracted  from  the  final  error  to  determine  the  actual 
error  in  a  particular  problem. 

2. 11.     REDUCTION  IN  ERRORS 

It  is  apparent  that  errors  could  be  reduced  if  the  size  of  the  limit 
cycle     is  reduced.     Two  approaches  will  be  taken  to  this  problem.    First, 
system  parameters  will  be  adjusted  in  an  attempt  to  minimize  the  limit 
cycle.     Second,   a  system  of  compensation  by  switching  logic  will  be 
investigated. 
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Figure  2-3.     Steady  State  Limit  Cycle:    T      =  1,  r  =  0. 
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Figure  2-5.     Steady  State  Limit  Cycle:    T     =  1,   r  =  0.2. 
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Figure  2-6.     Steady  State  Limit  Cycle:    T     =  1,  r  =  0.  3, 


Figure  2-7.     Steady  State  Limit  Cycle:    T     =  1,  r  =  0.  4. 
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Figure  2-8.     Steady  State  Limit  Cycle:    T     =  1,  r  =  0.  5. 
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Figure  2-9.     Steady  State  Limit  Cycle:    T     =  1,  r  =  0.6, 
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Figure  2-10.     Steady  State  .Limit  Cycle:    T     =  1,  r  =  0.7. 
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Fig.   2-11.     Steady  State  Limit  Cycle:    T     =  1,  r  =  0.  8, 
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Figure  2-12.     Steady  State  Limit  Cycle:    T     =  1,  r  =  0.  9, 
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CHAPTER  3 

MINIMIZING  LIMIT  CYCLE  SIZE  BY  PROPER  SELECTION 
OF  SYSTEM   PARAMETERS 

3.1,     GENERAL 

In  the  previous  chapter  it  was  shown  that  a  ternary  pulsed  pen- 
dulous accelerometer  may  operate  in  several  modes  of  oscillation  about 
a  bias  position.     The  bias,   or  average  pendulum  position,   produces  two 
distinct  errors.     First,    the  bias  itself  represents  an  error  in  velocity 
indication  in  that  a  certain  amount  of  acceleration  reaction  is  opposed 
by  viscous  friction  in  the  process  of  the  pendulum  moving  from  its  null 
position.     This  error  is  non-accumulative  and  its  effect  may  even  be 
negligible.     Second,    cross-coupling  error  is  introduced  by  the  bias 
angle.     For  small  angle  approximations,    the  cross-coupling  error  is 
the  time  integral  of  the  product  of  the  bias  angle  and  the  component  of 
acceleration  along  the  reference  axis,     Cross-coupling  error  is  ac- 
cumulative (i.  e.    it  may  grow  with  time).     It  may  be  minimized  by  main- 
taining the  bias  angle  small,   or  may  be  "computed  out"  if  the  bias  angle 
is  known. 

The  additional  error  associated  with  the  oscillation  of  the  pendu- 
lum about  its  bias  position,   is  oscillatory  and  therefore  non-accumulative. 
The  instantaneous  error  caused  by  this  oscillation  might  even  be  negli- 
gible.    However,   these  oscillations  if  large  may  cause  the  velocity  in- 
formation to  be  presented  to  the  computer  in  an  undesirable  form. 
Additional  errors  may  be  generated  by  the  computer  which  are  ac- 
cumulative.    The  nature  and  extent  of  these  errors  will  not  be  investi- 
gated here.     To  minimize  these  errors  it  is  desirable  to  maintain  the 
pendulum  in  its  lowest  mode  of  oscillation.     (A  limit  cycle  L.j  of  the 
lowest  mode  is  defined  to  be  that  limit  cycle   which  encloses  the  least 
area  of  the  non-dimensional  phase  plane.     Thus  for  a  non-dimensional 
acceleration  of  2/75   a  1:2:1:3  limit  cycle  is  of  lower  mode  than  a  2:5  or 
a  1:1:1:4  limit  cycle.  ) 

It  is  presupposed  that  certain  characteristics  are  specified  before 
initial  design  of  a  ternary  pulsed  accelerometer  begins.     Sample  time, 
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maximum  acceleration  and  physical  size  and  weight  limitations  may  be 
among  the  specifications.     The  designer  is  at  liberty  to  select  other 
parameters  in  proper  combination  to  insure  adequate  performance.    Such 
selection  must  be  compatible  with  physical  realizability  at  a  "reasonable" 
cost. 

It  has  been  previously  determined  that  dynamic  behavior  of  a 
ternary  system  is  determined  by  two  non-dimensional  quantities: 
sample  time,   T   ,   and  dead  zone,  $..     Knowledge  of  the  effect  of  these 
two  quantities  on  system  performance,  offers  the  designer  a  "crutch" 
in  design  analysis. 

3.  2.     SELECTION  OF  NON-DIMENSIONAL  SAMPLE  TIME 

Non-dimensional  sample  time  is  the  ratio  of  sample  time  to  the 
characteristic  time  of  the  viscous  damped  pendulum.     It  is  desirable  to 
select  system  parameters  such  that  limit  cycles  of  higher  order  are 
prevented.     Since  the  value  of  non-dimensional  sample  time  has  a 
definite  effect  upon  the  modes  of  oscillation  possible,   it  should  be 
selected  to  prevent  as  many  higher  order  limit  cycles  as  possible. 

Figures  C-2  through  C-6  show  the  freedom  in  position  of  a  simple 
limit  cycle  about  a  switching  line,   plotted  as  a  function  of  non-dimensional 
sample  time.     A  simple  limit  cycle  is  defined  to  be  a  limit  cycle  that  is 
complete  after  two  switches  between  torque  modes.     Appendix  C  points 
out  that  if  the  freedom  in  position  of  a  simple  limit  cycle  is  negative^ 
then  the  simple  limit  cycle  under  consideration  cannot  exist.     In  general, 
as  the  non-dimensional  sample  time,   T   ,  is  increased,  more  simple 
limit  cycles  are  prevented.     It  is  observed  that  the  2:2  limit  cycle  can- 
not be  prevented  by  adjusting  the  non-dimensional  sample  time  as  the 
freedom  in  position  of  the  limit  cycle  is  positive  for  all  values  of  sample 
time.     Likewise  the  1:N  limit  cycles,   where  N  is  any  integer,   cannot  be 
prevented.     The  1:N  limit  cycles  are  of  the  lowest  mode  and  are  desirable. 

A  value  of  five  or  greater  for  non-dimensional  sample  time  will 
prevent  all  non-desirable  simple  limit  cycles  with  the  single  exception 
of  the  2:2  limit  cycle.     It  is  realized  that  the  simple  limit  cycles  rep- 
resent only  a  few  of  the  infinity  of  limit  cycles  possible.     However,  one 
would  expect  that  the  dynamic  behavior  of  the  pendulum  would  not  change 
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greatly  for  a  small  change  in  acceleration.     For  this  reason,   it  is  felt 
that  the  simple  limit  cycles  are  representative  of  the  over-all  behavior 
of  the  pendulum. 

3.  3.     SELECTION  OF  NON-DIMENSIONAL  DEAD  ZONE 

The  non-dimensional  dead  zone  must  be  larger  than  the  non- 
dimensional  sample  time  in  order  to  prevent  torquing  in  a  direction  to 
aid  acceleration.     If  it  is  larger  than  necessary,   the  bias  angle  will  be 
larger  than  necessary  and  errors  associated  with  the  bias  angle  will  be 
larger.     It  would  appear  that  the  best  value  for  the  non-dimensional 
dead  zone  is  one  which  is  larger  than  the  non-dimensional  sample  time 
by  an  infinitesimal  amount. 

3.  4.     SELECTION  OF  DIMENSIONAL  PARAMETERS 

The  physical  dead  zone  that  is  inherent  to  the  error  detector  is 
a  limiting  factor  in  maintaining  the  null  error  small.     It  would  appear 
advisable  that  the  dead  zone  not  be  made  larger  artificially,   in  order 
to  implement  a  ternary  pulsed  system. 

The  key  to  the  proper  selection  of  dimensional  parameters  lies  in 
the  dimensional  to  non-dimensional  transformations  of  Table  2-1.     If  the 
physical  dead  zone  of  the  error  detector  is  taken  to  be  the  value  that  will 
transform  into  a  non-dimensional  dead  zone  that  is  equal  to  the  non- 
dimensional  sample  time,   a  system  of  equations  is  introduced  which 
may  be  used  in  the  selection  of  dimensional  parameters.     The  fact  that 
the  non-dimensional  dead  zone  must  be  larger  than  the  non-dimensional 
sample  time  by  a  small  amount,   can  be  taken  care  of  by  final  adjust- 
ment to  the  threshold  detector. 

Suppose  that  A*    is  the  dead  zone  expected  from  the  error  de- 
tector.    Then  by  the  transformations  of  Table  2-1: 

AAC2 

<K=  — (3.1) 

^  MJ 

and 

Ct 

T     =  — -   .  (3.2) 

s 
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By  desire  T  and  </).  are  set  equal 


Ts=*A.  (3.3) 

The  applied  torque,   M,  must  be  sufficient  to  overcome  the  reaction 
of  the  maximum  expected  acceleration.     Therefore: 

M  =  a  P.  (3.4) 

max 

If  a  ,  AA,   and  t     are  specified,   six  parameters  remain  to  be 

max'      A  s  M  *  r 

fixed  by  the  designer  to  satisfy  the  above  four  equations.     The  designer 
may  select  a  value  for  T     to  give  satisfactory  dynamic  performance, 
leaving  five  parameters  to  be  adjusted  to  satisfy  the  four  equations. 

3.  5.     EXAMPLE 

Suppose  it  is  desired  to  design  a  ternary  pulsed  pendulous  ac- 

4  2 

celerometer  with  a  maximum  acceleration  capability  of  10     cm/ sec     and 

with  a  sample  time  of  500  microseconds.     The  velocity  quantum  is  the 
product  of  these  two  quantities  or  5  cm/  sec.     A  value  of  five  is  selected 
for  non-dimensional  sample  time  so  that  all  undesirable  simple  limit 
cycles  are  prevented  with  the  exception  of  the  2:2  limit  cycle.     Suppose 
also  that  the  dead  zone  expected  from  the  error  detector  is  five  micro- 
radians.     Substituting  these  values  into  Equations    (3.1)  through  (3.4): 

—  =  100  sec.  (3. 5) 

M 


—   =  104/sec  (3.6) 


M  =  104cm/sec2.  (3.7) 


If  P  is  arbitrarily  selected  to  be  one  gm-cm.,   Equations  (3.  5) 
through  (3.7)  can  then  be  solved  for  the  remaining  unknowns: 
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4 
M  =  10     dyne-cm 

2 

J  =  100  gm-cm 

C  =  10     dyne-cm-sec. 

Figures  3-1  through  3-11  show  the  steady  state  phase  plane  response 
of  such  a  system  for  various  values  of  input  acceleration.     Figures  3-12 
through  3-15  show  the  null  error  versus  time  response  for  the  phase 
plane  responses  of  Figures  3-1  through  3-11. 

The  phase  plane  plots  show  that  the  pendulum  is  in  its  lowest  mode 
of  oscillation  except  when  the  acceleration  is  in  the  vicinity  of  one  half 
the  maximum.     This  could  be  expected  in  that  the  2:2  limit  cycle  has  not 
been  eliminated.     If  the  2:2  limit  cycle  does  not  produce   excessive  error 
in  the  computer,   or  if  the  accelerometer  is  not  expected  to  operate  in  the 
vicinity  of  one  half  the  maximum,   then  the  system  is  satisfactory  in  that 
sufficient  higher  order  limit  cycles  have  been  eliminated. 

3.6.     LIMITATIONS  TO  THE  METHOD 

The  authors  admit  that  the  previous  example  has  been  somewhat 
"rigged"  so  that  realistic  answers  could  be  obtained  for  dimensional 
parameters.     The  specifications  could  be  such  that  Equations  (3.1) 
through  (3.  4)  produce  dimensional  quantities  that  are  not  physically 
realizable  if  a  value  of  five  is  selected  for  T   .     One  may  have  to  settle 
for  a  value  much  less  than  five  for  T   ,   in  which  case  another  method 
will  have  to  be  employed  to  maintain  the  limit  cycles  in  their  lowest 
mode.     Such  methods  are  the  subject  of  the  next  chapter. 
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Figure  3-13.     Null  Error:    T     =  5,   r  =  0.  4,    0,45,   0.5. 
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Figure  3-14.     Null  Error:      T     =  5,  r  =  0.55,   0.6,   0.7, 
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Figure  3-15.     Null  Error:    T     =  5,   r  =  0.  8,   0.9. 
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CHAPTER  4 

COMPENSATION  OF  TERNARY  MODE  PULSE  TORQUED 
PENDULOUS  INTEGRATING  ACCELEROMETERS 

4.1.     INTRODUCTION 

Torquing  only  in  the  direction  to  oppose  the  acceleration  reaction 
torque  and  torque  pulsing  in  a  routine  of  the  lowest  order  are  the  en- 
visioned benefits  of  ternary  mode  pulse  torquing.     By  routine  of  the 
lowest  order  is  meant  torquing  with  a  routine  which  keeps  the  area  en- 
closed by  the  limit  cycle  as  small  as  possible.     When  the  input  accel- 
eration-maximum acceleration  ratio  r,   is  less  than  one-half,   this  will 
be  done  if  the  restoring  torque  pulse  is  the  shortest  time  possible,  which 
is  T   ,   the  non-dimensional  sampling  period.     For  values  of  r  greater 
than  one-half,   the  non-torquing  period  should  have  the  shortest  possible 
time  duration,   again  T   . 

Any  rational  value  of  r  can  be  represented  by  a  fraction  with 
integer  numerator  and  denominator. 

r  =  -*-,  (4.1) 

K+  L 

where  K  is   the  total  number  of  torquing  sample  periods  per  counting 
cycle,   and  L  is    the  total  number  non-torquing  sample  periods  per 
counting  cycle.     A  counting  cycle  is  of  K  +  L  sample  periods  duration 
and  is  the  smallest  number  of  sample  periods  which  accurately  rep- 
resents r.     Both  K  and  L  can  be  the  number  of  sampling  periods  in 
more  than  one  torque  or  non- torque  periods.     Or 

K  =  KL  +  K2  +  K3  +    (4.  2} 

where  K,,  K2,  K~,   ••  •  *   are  the  number  of  sampling  periods  in  individual 
torque  periods,   and  i 

L  =  Lx  +  L2  +  L3  + (4.  3) 

where  L,,   L?,   L„,   •  •  •  •    are  the  number  of  sampling  periods  in  individual 
non-torque  periods.     When  r  is  less  than  one-half  the  lowest  order  counting 
cycle  can  be  produced  if  K,,   K2,   Kg,    •  •  •  •    equal  one.     Likewise,   when  r  is 
greater  than  one-half,   the  lowest  order  counting  cycle  can  be  produced  if 
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L^,   L2,   Lo,   ••••    equal  one. 

The  above  criteria  will  be  the  goal  in  compensating  pulse  torqued 
pendulous  accelerometers.     In  achieving  this  goal  errors  other  than 
the  error  due  to  higher  order  counting  cycles  will  be  introduced.     These 
new  errors  should  be  smaller  than  the  errors  produced  by  the  larger 
limit  or  counting  cycles,  or  they  must  be  of  a  nature  such  that  the 
total  velocity  stored  in  the  computer  can  be  corrected  for  the  errors 
introduced  by  the  compensation. 

The  larger  order  limit  cycles  found  in  sample  data  bang-bang 
servos  are  the  result  of  lag  introduced  in  the  system  by  discrete  sam- 
pling rather  than  continuous  sampling.     Thus,   lead  must  be  introduced 
into  the  error  circuit.     In  nonlinear  servo  theory,  on  the  error-error 
rate  phase-plane,  a  servo  with  lead  in  the  error  circuit  and  a  dead 
zone  has  its  switching  lines  represented  in  the  manner  shown  in  Fig- 
ure 4-1'.     In  the  case  of  sample  data  bang-bang  servos,   switching  from 
one  torque  state  to  another  occurs  at  the  first  sampling  instant  after 
crossing  the  switch  line. 

In  servo  mechanisms  in  which  the  error  magnitude  and  sense  is 
represented  in  an  analog  manner  by  a  voltage,   lead  compensation  can 
be  inserted  by  active  or  passive  electronic  networks.     In  pulse  torqued 
integrating  accelerometers  null  error  is  analoged  by  the  voltage  output 
of  the  microsyn.     Representative  sampling  frequencies  are  one  to  five 
kilocycles  (1000  to  5000  c.  p.  s.  ).     Maximum  frequencies  at  which  micro- 
syns  have  been  found  to  be  excited  are  also  one  to  five  kilocycles.     Phase 
distortion  increases  as  the  exciting  frequency  is  increased.     The  error 
analog  is  of  poor  quality.     In  existing  binary  pulse  torqued  pendulous 
accelerometers,    the  microsyn  null  position  error  detector  is  used 
only  to  determine  the  sense  of  the  error.     Hence,  because  of  poor  error 
analog  quality,   compensation  of  the  error  analog  signal  is  not  practical. 

What  is  available  in  the  system  is  the  error  sense  sampled  at 
equal  periods.     This  suggests  compensation  using  logic  techniques. 
To  simulate  the  lead  switch  lines  shown  in  Figure  4-1,   error  sense, 
magnitude,   and  rate  must  be  obtainable.     The  best  magnitude  resolution 
which  can  be  expected  from  a  microsyn  type  position  error  detector  is 
about  the  width  of  the  microsyn  dead  zone.     Error  rate  is  not  readily 
available.     The  most  which  can  be  obtained  readily  from  microsyn  null 
error  detectors  at  the  desired  sampling  frequency  is  pendulum  position 

-67- 


+  Error  Rate 


Positive 
Torque 
Zone 


Dead 
Zone 


CUf 


SWITCH  LINE^Z-? 


Negative 
Torque 
Zone 


+  Error 


■■^-SWITCH  LINE 


Figure  4-1.     Phase  Plane  Representation  of  a  Second  Order 
Position  Servo  With  Dead  Zone  and  Lead. 
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in  zones.     These  zones  are  the  dead  zone  and  magnitude  zones  of 
detectable  sense. 

The  above  zones  can  be  determined  and  controlled  in  width  with 
limited  ability  by  passing  the  microsyn  output  signal  through  a  phase 
detector  to  determine  sense.     Then  the  signal  can  be  applied  to  biased 
level  detectors  to  set  zone  magnitude. 

Figure  4-2  shows  the  phase  plane  divided  into  five  zones. 
Following  are  logic  techniques  in  using  these  zones  to  insert  lead 
compensation  into  the  system. 

4.  2.     ALTERNATING  PULSE  ZONE  COMPENSATION 

Lead  compensation  can  be  interpreted  as  a  means  of  anticipating 
required  actions  and  carrying  them  out  sooner  than  a  non-compensated 
system  would.     Figure  4-1  shows  that  the  magnitude  of  the  lead  effect 
should  be  proportional  to  the  magnitude  of  the  velocity.     This  is  not  pos- 
sible with  available  error  detectors  at  the  information  rate  required. 
Where  our  system  does  not  allow  a  pure  lead  effect,   the  lead  producing 
zone  can  produce  a  quasi  lead  effect. 

A  zone  of  sampling  periods  of  alternating  torquing  and  not  torquing 
periods  will  anticipate  both  the  change  from  the  non- torquing  dead  zone 
to  the  continuous  torquing  zone  and  vice  versa.     Considering  positive 
input  accelerations,   it  is  obvious  for  values  of  input  acceleration  - 
maximum  acceleration  ratio,   r,   greater  than  0.  5  that  more  torquing 
periods  are  required  in  a  counting  cycle  than  non-torquing  periods. 
Hence,  </>  will  increase  until  switching  is  about  the  boundary  between 
the  continuous  and  alternating  zones,   the  -&•  line.     Similarly  for  values 
of  r  less  than  0.  5  the  switching  will  be  about  the  boundary  between  the 
alternating  zone  and  the  dead  zone,   the  -<£pvline. 

Previously  the  non-dimensional  equations  of  motion  of  the  pendulum 
were  determined  to  be 

<j>  =  (6  -  r)(T  +  e~  '  -  1)  +(1  -  e'T)i     +  <f>  (4.4) 


i  =  (6  -  r)(l  -  e    T)  +iQe"T  (4.5) 
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Except  for  the  situation  in  which  a  very  quick  change  in  the  accelera- 
tor has  reversed  the  sense  of  the  acceleration  such  that  the  accelera- 
tion reaction  torque  and  the  applied  error  reduction  torque  are  in  the 
same  direction,   the  magnitude  of  maximum  $   in  the  alternating  torquing 
zone  is  seen  from  equation  (4.  2)  to  be  (6  -  r).     One-half  is  the  maxi- 
mum value  of  r  when  switching  from  the  </>n  line.     Since  when  moving 
from  the  dead  zone  into  the  alternating  torque  zone  6  equals  zero,  the 
maximum  </>  associated  with  switching  about  the  -$n  line  in  the  alter- 
nating torque  zone  is  0.  5.     Similarly,   since  0.  5  is  the  minimum  value 
of  r  and  6=1  when  switching  about  the  -<j>T    line  and  moving  from  the 
continuous  torque  zone  into  the  alternating  torque  zone,   the  maximum 
magnitude  of  </>  associated  with  switching  about  the  -</>,    line  in  the 
alternating  torque  zone  is  0.  5. 

If  <j>  is  greater  than  0.  5  in  the  alternating  torque  zone  the 
switching  locus  is  moving  from  one  boundary  to  the  other  because  of 
a  change  in  r  from  larger  than  0.  5  to  smaller  than  0.  5,  or  vice  versa. 
In  this  case  switching  in  the  alternating  zone  is  not  necessary.     If  </>  is 
less  than  0.  5  switching  is  mandatory.     This  will  be  assured  if  the 
width  of  the  alternating  zone  is  made  at  least  equal  to  the  maximum  $ 
desired  to  be  "caught"  in  the  alternating  zone  times  the  non-dimensional 
sampling  period  T   .     The  minimum  value  of  </>_    -  <£_  acceptable  is 

S  L  D 

0.  5  T   . 
s 

The  width  of  the  dead  zone  must  be  adjusted  so  that,  for  very 
small  values  of  r,   the  pendulum  does  not  pass  into  the  alternating  zone 
where  the  magnitude  of  (6  -  r)  is  greater  than  1.     Refer  to  the  conditions 
shown  in  Figure  4-3  for  the  limiting  case  of  small  r,   that  is,  r  approxi- 
mately  equal  to  zero.     </>     equals  -</>n>  (j>     approximately  equals  zero. 
</>2  is  less  than  +</>r).     </>2  is  the  point  of  furthest  penetration  into  the 
positive  </>  plane  before  the  pendulum  is  driven  back  to  the  negative  </> 
plane.     By  equations  (4.  4)  and  (4.  5),  with  a  torque  pulse  of  magnitude 
(6  -  r)  approximately  equal  1,   the  position  one  sampling  period  later  is 

-T 
^  =  Tg+e       S-l-(/>D  (4.6) 

-T 

6,   --   1  -  e       S  (4.7) 
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Figure  4-3.     Phease  Plane  Condition  for  Determining  Minimum  Width 
of  the  Dead  Zone  in  a  Compensated  Ternary  Pulsed  In- 
tegrating Pendulous  Accelerometer. 
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When  the  pendulum  is  coasting,    that  is,    6  equals  zero  and  r  equals 
zero,  and  the  phase  plane  path  has  not  crossed  the  zero  velocity  axis, 
the  motion  can  be  described  by 

*+i=*0+i0.  (4.8) 

Hence,   when  </>  equals  zero  at  </»2 

</>2  =  <j>l+  </>1  (4.9) 

-T  -T 

<j)2  =  Tg  +  e    S  -  1  -  </>D  +1  -  e      s  (4.10) 

^2  =  Ts  "  *D*  (4-U) 

To  keep  §  in  the  dead  zone, 

(+<j>D)  -  (-4>D)  >  4>2  -  (-</>D),  (4.12) 

(+<j>D)  -  (-*D)  >  Tg,  (4.13) 

or,  the  non-dimensional  magnitude  of  the  minimum  width  of  the  dead 
zone  is  equal  to  the  magnitude  of  the  non-dimensional  sampling 

period,   T   . 

s 

An  elemental  switching  logic  is  explained  and  designed  in 
section  D.  2  of  Appendix  D.     A  symbolic  design  is  shown  in  Figure  D-4. 

Figures  (4-4),  (4-5),  and  (4-6)  show  the  steady  state  counting 
cycles  for  the  alternating  pulse  zone  compensation  for  values  of  r  in 
one-tenth  increments  from  0.1  to  0.9.     These  values  were  chosen 
because  the  counting  cycles  are  relatively  simple  yet  they  exhibit  the 
problems  found  in  the  technique.     Where  practical,   two  counting 
cycles  are  shown  and  they  indicate  the  boundaries  of  the  area  in  the 
phase  plane  in  which  the  limit  cycles  can  exist. 

For  the  cases  of  r  equal  0.  3  and  0.  7  only  one  limit  cycle  is 
shown.     Inspection  of  the  cycles  show  they  have  a  very  limited  freedom 
of  position,  </>».     The  trajectory  of  the  error  on  the  phase  plane  will 
adjust  itself  and  hunt  for  the  limit  cycle.     In  Figure  (4-4)  the  area  in 
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which  the  trajectory  for  r  equal  0.  3  hunts,   after  large  transient  effects 
have  died  out,   is  shown  with  dashed  lines.     One  can  see  that  the  over- 
all counting  error  during  the  hunting  period  is  about  that  of  the  steady 
state  limit  cycle. 

For  values  of  r  less  than  0.  5  the  counting  cycles  are  found  to  be 
in  the  lowest  order.     For  values  of  r  greater  than  0.  5  requiring  a  two 
period  torque  pulse  the  counting  cycles  are  not  in  their  lowest  order, 
r  equal  0.6  requires  a  torque  2  -  not  torque  1,   torque  3  -  not  torque  1 
counting  cycle  for  lowest  order  instead  of  3  -  1,   1-1,   1-1,  1-1  counting 
cycle,     r  equal  0.  7  requires  3-1,   2-1,   2-1  counting  cycle  instead  of 
5-1,   1-1,   1-1  counting  cycle. 

In  order  to  have  a  two  period  pulse  in  the  alternating  pulse  zone 
compensation,   the  pulse  must  start  in  the  continuous  positive  torque 
zone  and  move  into  lead  compensation  zone  before  the  end  of  one  sam- 
pling period  in  time.     For  different  values  of  r,   meeting  the  above  con- 
dition depends  on  the  value  of  T   ,   the  non-dimensional  sampling  time, 

to 

or  the  ratio  of  the  sampling  period  to  the  time  constant  of  the  pendulum. 

Figures  4-7  through  4-12  show  the  steady  state  counting  cycles 
for  alternating  pulse  zone  compensation  with  T     equal  5.  0.     Those 

to 

values  of  r  which  were  in  their  lowest  order  at  T     equal  1.  0  have  a 

s 

larger  freedom  of  position  at  T     equal  5.  0.     For  r  equal  0.6, 

to 

Figure  4-10,   the  lowest  order  counting  cycle  can  now  be  assumed  with 
some  freedom  of  position,  $».     For  r  equal  0.7,   Figure  4-9,   the  lowest 
order  counting  cycle  is  possible,  but  the  system  is  required  to  switch 
at  almost  a  point  on  the  phase  plane.     The  system  will  "hunt"  for  this 
point  in  the  3-1,   1-1,    3-1  order  after  large  transient  effects  have 
died  out. 

Figure  4-11,  for  r  equal  0.  88  and  T     equal  5.  0,   shows  the 

to 

counting  cycle  is  in  its  lowest  order;     8  -  1,   7  -  1,   7  -  1.     At  T     equal  1.  0 

to 

the  <f>„  is  approximately  one -fifth  that  shown  in  Figure  4-12,  again  re- 
quiring more  "hunting"  to  get  in  a  lowest  order  counting  cycle.     For 
larger  values  of  r  which  do  not  require  a  two  sample  period  pulse  in 
their  lowest  counting  order,   the  lowest  order  is  possible  for  small 
values  of  T   ,  but  increasing  T     reduces  the  amount  of  "hunting"  re- 

to  to 

quired  to  assume  the  lowest  order  counting  cycle. 
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4.  3.     INVERSE  ORDER  ALTERNATING  PULSE  ZONE  COMPENSATION 

In  the  alternating  pulse  zone  compensation  enough  lead  effect 
was  introduced  as  the  position  entered  the  lead  compensation  zone  from 
the  dead  zone.     However,   on  entering  the  lead  compensation  zone  from 
the  continuous  torque  zone  the  system  continues  torquing  one  sample 
period  after  the  first  sample  in  the  lead  compensation  zone  before  it 
switched  to  non- torquing.     This  removed  most  of  the  lead  effect.    By 
changing  the  switching  logic  so  that  when  a  sample  is  in  the  lead  com- 
pensation zone  the  action  during  the  period  following  the  sampling 
instant  is  opposite  the  action  carried  out  the  previous  period,    that  is 
non-torquing  if  torqued  previous  period  and  vice-versa,   an  inverse 
order  alternating  pulse  zone  is  set  up.     This  will  produce  the  desired 
lead  effect  in  both  the  increasing  and  decreasing  <f>  direction. 

For  the  same  reasons  as  in  the  alternating  pulse  zone  compen- 
sation,   the  width  of  the  dead  zone  must  be  greater  than  T     and  the  width 
of  the  lead  compensation  zone  greater  than  0.5  T   , 

The  switching  logic  for  the  inverse  order  alternating  pulse  zone 
compensation  is  developed  in  section  D-3  of  Appendix  D  and  shown  in 
the  diagram  of  Figure  D-5. 

The  lead  effect  when  switching  about  the  (fry.  lines  was  not 
changed.     The  phase  plane  results  for  values  of  r  from  0  to  0.  5  are 
the  same  as  for  the  alternating  pulse  zone  compensation  shown  in 
Figures  4-4  and  4-5.     The  phase  plane  results  for  T     equal  1  and  r 
greater  than  0.  5  are  shown  in  Figure  4-12  and  4-13.     The  resulting 
lowest  order  counting  cycles  are  shown.     Also  Figure  4-13  shows  the 
expected  result,   that  if  T     is  increased  the  counting  cycles  remain 

fa 

in  their  lowest  order. 

If  the  compensation  problem  is  viewed  from  the  aspect  of  not 
allowing  higher  order  limit  cycles  by  parameter  variation  as  de- 
veloped in  Chapter  3,   the  same  results  are  obtained  as  in  the  inverse 
order  alternating  pulse  zone  compensation  except  as  to  the  width  of 
the  compensation  zone. 

Referring  to  Figure  C-l,   if  the  portion  of  the  switching  line  for 
positive  angular  rate,   $  ,   is  moved  to  the  left  of  its  present  position, 
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the  N  =  K  -  1  sample  will  be  to  the  right  of  the  switch  line  and  the 
system  will  stop  torqui.ng.     Similarly  if  the  portion  of  the  switching 
line  for  negative  angular  rate  is  moved  to  the  right  of  its  present 
position,   the  N  =  K  +  L  -  1  sample  will  be  to  the  left  of  the  switch 
line  and  the  system  will  start  to  torque.     The  action  described  has 
the  effect  cf  making  the  freedom  in  positon,  <f>~,  negative  for  higher 
order  limit  cycles. 

For  the  lowest  order  limit  cycles,   and  r  greater  than  0.  5,   the 
system  may  non- torque  only  one  sample  period,   T   ,   at  a  time.     The 
positive  </>   portion  of  the  switching  line  can  best  be  instrumented  by 
the  positive  portion  of  the  </>,    switching  line  of  the  inverse  order 
alternating  compensation  zone.     By  similar  reasoning,   the  negative 
$  portion  of  the  switching  line  can  best  be  instrumented  by  negative 
portion  of  the  </>n  switching  line  of  the  inverse  order  alternating  com- 
pensation zone.     Subsequent  samples  still  in  the  compensation  zone 
will  cause  alternating  torquing. 

From  Figure  C-2  we  see  if  the  switching  line  separation  is  made 
greater  than  0.  363  for  T     equal  1.  0  the  2:2  limit  cycle  is  eliminated. 
Figures  C-3  through  C-6  show  if  <£,,  eliminates  the  2:2  limit  cycle, 
all  higher  limit  cycles  are  eliminated.      |</>T    -  <£_J   effects  </>,..     Inspec- 
tion of  Figure  C-2  shows  that  as  T     increases  the  required  $,.  does 
not  increase  linearly.     In  fact  it  does  not  exceed  1.  0. 

The  logic  required  is  the  logic  for  the  inverse  order  alter- 
nating torque  zone  compensation.     Viewing  the  compensation  problem 
from  the  parameter  variation  viewpoint  resulted  in  the  same  com- 
pensation logic  but  with  the  minimum  width  of  the  dead  zone  decreased 

from  0.  5  T     to  0.  363  T    . 

s  s 

The  result  of  decreasing  the  width  of  the  lead  compensation  zone 
is  to  decrease  the  error  bias  and  cross  coupling  effects.     But  for  very 
large  transient  situations  the  time  required  to  reach  a  steady  state 
limit  cycle  will  be  longer. 

The  phase  plane  trajectories  for  T     equal  one  are  similar  to 
Figures  4-12  and  4-13  except  that  the  -</>T    line  is  at  -0.  363  instead  of 
-0.  5  and  thus  the  error  bias  is  reduced  by  0.137. 
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4.  4.      LOGICALLY  FORCED  LOWEST  ORDER  COMPENSATION. 

In  this  compensation,   switching  logic  is  devised  to  allow  the  sys- 
tem to  torque  only  one  period  at  a  time  if  r  is  less  than  0.  5  and  to  not 
torque  only  one  period  at  a  time  if  r  is  greater  than  0.  5.     Also  limit 
cycle  switching  will  be  limited  to  be  about  the  </>„  line  in  order  to  re- 
duce cross  coupling  errors.     If  r  is  greater  than  0.  5  and  the  system 
only  torques  one  period  at  a  time  the  limit  cycle  position  will  move 
toward  greater  error.     If  r  is  less  than  0.  5  and  the  system  will  not 
torque  only  one  period  at  a  time,   the  limit  cycle  position  will  move 
to  smaller  error. 

Based  on  the  described  movement,   when  it  is  desired  to  switch 
from  logic  for  r  less  than  0.  5   to  logic  for  r  greater  than  0.  5  the 
position  will  be  moving  toward  greater  error.     If  a  zone  or  boundary 
line  called  +$T    or  -$,    as  appropriate,  is  set  up  beyond  the  greatest 
error  of  counting  cycles  for  r  less  than  0.  5,   when  r  is  greater  than 
0.  5  the  position  will  increase  until  it  passes  this  line.     When  a  sam- 
ple with  <j)  greater  than  $,    is  obtained  the  logic  can  be  changed  to 
torque  as  required  when  out    of  the  dead  zone,   but  not  torque  less 
than  one  sampling  period,   T   ,   at  a  time.     The  system  will  then  torque 
until  the  position  returns  to  the  dead  zone  and  then  switching  will  be 
about  the  <f>j^  line.     This  logic  system  requires  the  position  to  return 
out  of  the  dead  zone  in  only  one  period  of  not  torquing.     To  insure  this 
for  all  combinations  of  dead  zone  penetration  and  all  values  of  r 

greater  than  0.  5  it  is  found  that  T     must  be  at  least  5. 

°  s 

The  dead  zone  width  requirements  are  the  same  as  before, 
(+</>)  -  (-</>)  greater  than  T   . 

LJ  S 

Then  with  T     equal  5  or  greater,   r  will  be  less  than  0.  5  when 
two  consecutive  samples  are  found  to  be  in  the  dead  zone.     This  can 
be  used  to  change  the  switching  logic  when  r  becomes  less  than  0.  5. 
For  consecutive  samples  in  the  zone  between  </>,    and  <£„,   the  system 
will  torque  alternating  periods. 

The  logic  for  this  compensation  is  developed  in  section  D-4  and 
shown  in  Figure  D-6. 

Figures  4-14  through  4-18  show  the  counting  cycle  for  the 
logically  forced  lowest  order  compensation  to  be  in  the  lowest  order. 
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Figure  4.-14.     Steady  State  Limit  Cycles  for  Logically  Forced 
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In  fact  for  values  of  r  other  than  in  the  0.  5  vicinity  the  counting  cycles 
are  exactly  the  same  as  in  the  non-compensated  system  with  T     equal 
5.  0.     Only  in  the  vicinity  of  0.  5  is  the  counting  cycle  improved. 

4.  5.     COMPENSATION  RESULTS 

On  each  of  the  figures  showing  limit  cycles  the  zero  error  line 
is  shown.     This  allows  one  at  a  glance  to  get  a  relative  estimate  of 
the  error  size.     In  Figure  4-19  through  4-24  the  values  of  the  error 
during  the  limit  cycles  are  shown.     The  errors  were  computed  using 
the  previously  developed 

E  =  <j>  +  <j>   .  (4.14) 

E  is  the  non-dimensional  error  expressed  in  counts  of  velocity  quanta. 
The  largest  error  allowed  by  the  freedom  of  position  is  shown.     The 
following  conclusion  can  be  drawn  from  the  error  plots  and  the  phase 
plane  trajectories. 

In  the  alternating  pulse  zone  compensation  T     must  be  in- 
creased  to  5  to  insure  the  counting  cycles  are  in  their  lowest  order. 
The  only  performance  improvement  is  in  the  vicinity  of  r  equal  0.  5. 
The  hardware  restrictions  placed  on  the  system  by  the  uncompensated 
system  due  to  a  large  T     has  not  been  removed.     In  addition,  for  values 
of  r  greater  than  0.  5  the  error  bias  has  doubled.     The  more  compli- 
cated error  detection  and  switching  logic  is  a  relatively  high  price  to 
pay  for  doubtful  performance  improvement. 

The  logically  forced  lowest  order  compensation  also  requires 
T     greater  than  5.     It  has  the  same  error  bias  as  the  uncompensated 
system.     Its  counting  cycles  are  the  same  as  those  of  the  uncom- 
pensated system    with  T     equal  5  except  when  riis  in  the  vicinity  of 
0.  5.     Theologically  forced  lowest  order  compensation  does  produce 
the  lowest  order  counting  cycles  for  all  values  of  r. 

The  inverse  order  alternating  pulse  zone  compensation  pro- 
duces the  lowest  order  counting  cycle  for  all  values  of  r  with  T     equal 
one  or  greater.     For  r  equal  0.  5  and  greater  the  error  bias  is  greater 
than  in  an  uncompensated  system  due  to  the  switching  about  the  </>. 
line.     The  switching  logic  shown  in  Appendix  D  for  the  inverse 
alternating  pulse  zone  compensation  is  relatively  simple. 
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Figure  4-19.     Error  Cycles  for  Compensated  Systems  With  T     =1.0 
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Figure  4-20.     Error  Cycles  for  Compensated  Systems  With  T     =1.0. 


■95- 


P 

O 

o 

< 

P 

« 
o 

H 
P 


-0.5 


-1.0 


-1.5 


-2.0 


t  ■  10  t 


r  =  0.  8 


=  10  t 


r  =  0.9 


Inverse  Alternating  |$      -  </>_|    =  0.  36  3  T 

Li  L)  S 

—  ^_ —  Inverse  Alternating  |<£_    -  <£_  |=0.5  T 

__  . .  Alternating 
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Figure  4-23.     Error  in  Compensated  Systems  With  T     =5.0, 


-98- 


0 
-0.5 
-1.0 
-1.5 


'\ 


V 


v 


A 


V 


s 


_S^ 


y\ 


10  t 


^\ 


H 


r  =  0.7 


O 

o 

< 

< 

O 

« 

w 

J 

13 


0.5 


■1.0 


-1.5 


■*\ 


41 


-^V 


10  t 


n; 


r  =  0.8 


=  10  t 


-2.  0 


r  =  0.  9 


Figure  4-24.     Error  in  Compensated    Systems 


With  T     =5.0. 
s 


-99- 


The  velocity  quanta  size,   the  maximum  acceleration  to  be  meas- 
ured,  and  the  computer  clock  frequency  usually  are  fixed  by  the  appli- 
cation of  the  pulsed  pendulous  accelerometer.     The  dead  zone  of  the 
error  detector  is  fixed  in  minimum  size.     These  fixed  parameters, 
with  the  additional  requirement  that  the  width  of  the  dead  zone  be 
greater  than  one  T   ,  place  stringent  demands  on  the  ratio  of  the  pen- 
dulum  moment  of  inertia  J  to  the  damping  constant  C.     Smaller  mini- 
mum required  T     eases  the  restriction  on  the  J  to  C  ratio. 

With  the  non-dimensional  width  of  the  dead  zone  equal  to  one  T   , 
the  error  bias  due  to  switching  about  the  </>.-.  line  is  one-half  a  velocity 
quanta  and  that  due  to  switching  about  the  (/>,    line  is  one  velocity 
quanta.     Sections  D.  5  and  D.  6  of  Appendix  D  develop  logic  to  correct 
the  velocity  count  for  the  error  bias  due  to  switching  about  particular 
switching  lines.     This  correction  makes  it  practical  to  make  the  dead 
zone  more  than  one  T     wide  in  non-dimensional  magnitude.     This  will 
ease  the  restrictions  put  on  the  J  to  C  ratio.     Correcting  the  velocity 
count  for  the  error  bias  does  not  decrease  the  physical  angle  the  pen- 
dulum makes  with  its  null  position.     It  will  not  decrease  cross  coupling 
effects. 

With  proper  error  bias  correction  and  using  the  inverse  order 
alternating  pulse  zone  compensation,  the  counting  error  can  be  kept 
within  one  quanta. 

In  describing  the  limitations  of  the  error  detector  the  conserva- 
tive statement  was  made  that  the  best  magnitude  resolution  to  be  ex- 
pected was  about  equal  to  the  width  of  the  dead  zone.     This  assumption 
was  based  on  the  high  noise  to  signal  ratio  due  to  a  relatively  low  error 
information  rate.     If  quadrature  effects  about  the  microsyn  null  posi- 
tion are  the  major  factor  in  creating  the  error  detector'  s  dead  zone 
size,  better  resolution  away  from  the  dead  zone  can  be  expected. 

Better  resolution  away  from  the  dead  zone  will  allow  the  non- 
dimensional  dead  zone  to  be  several  T     wide,  but  the  distance  be- 

s 

tween  </>,    and  </>      still  only  needs  to  be  0.  5  T   .    </>,    can  then  be  moved 
in  towards  the  dead  zone  decreasing  cross  coupling  effects.     This  will 
require  error  count  correction  because  the  dead  zone  will  become 
the  several  counts  wide.     Also,  because  the  dead  zone  is  several 
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counts  wide  the  random  initial  error  increases.     This  could  be  pre- 
vented by  electrically  restraining  the  pendulum  until  the  count  is 
initiated.     This  would  require  a  special  counting  error  bias  initial 
correction 

4.6.     COMPENSATION  SUMMARY 

Three  logic  compensation  schemes  were  devised.     The  alter- 
nating pulse  zone  compensation  with  a  small  amount  of  switching  logic 
traded  larger  error  bias  for  keeping  the  counting  cycle  in  its  lowest 
mode  for  all  values  of  r.     It  required  T     equal  to  five  or  more  placing 
design  restrictions  on  the  hardware.     The  logically  forced  lowest 
order  did  not  increase  the  error  bias  in  making  all  counting  cycles 
the  lowest  order,  but  it  again  requires  T     equal  five  or  more  and  uses 
more  complicated  switching  logic.     The  inverse  order  alternating  pulse 
zone  compensation  puts  the  counting  cycles  in  their  lowest  order  and 
only  requires  T     equal  to  or  greater  than  one.     The  switching  logic  is 
simple.     For  values  of  r  greater  than  0.  5  the  error  bias  is  increased. 
Velocity  count  correction  schemes  were  devised  to  correct  the  counted 
velocity  for  the  existing  error  bias.     The  inverse  order  alternating 
pulse  zone  compensation  is  the  most  practical  because  it  does  not 
place  stringent  moment  of  inertia  -  damping  coefficient  ratio  re- 
strictions on  the  pendulum  unit. 
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CHAPTER  5 
CONCLUSIONS  AND  RECOMMENDATIONS 

5.1.  CONCLUSIONS 

The  theoretical  feasibility  of  ternary  pulse  torquing  a  pendulous 
accelerometer  has  been  shown.     The  inherent  greater  stability  of  a 
ternary  system  as  compared  to  a  binary  system  results  in  smaller 
limit  cycles  associated  with  smaller  velocity  quanta  counting  error. 

By  adjusting  the  ternary  pulsed  accelerometer  system1  s  com- 
ponents so  as  to  make  the  non-dimensional  sample  time,   T  ,   5.  0  or 
greater,   the  lowest  order  limit  cycles  can  be  assured.     When  practical 
considerations  do  not  allow  adequate  adjustment,   switching  logic  com- 
pensation can  be  used  to  insure  the  lowest  order  limit  cycles  with  T 
equal  or  greater  than  1.  0.     In  compensating  the  ternary  mode  pulsed 
pendulous  accelerometer  increased  error  bias  is  traded  for  smaller 
limit  cycles.     If  the  error  bias  is  too  large,   the  velocity  quanta  count 
in  the  computer  can  be  corrected  for  the  error  bias.     Another  result 
of  the  compensation  is  increased  cross-coupling  effects.     These  can- 
not be  as  easily  corrected  in  the  velocity  count  and  thus  must  remain 
small. 

Ternary  mode  of  pulse  torquing  pendulous  integrating  accel- 
erometers  allows  increased  accelerometer  performance  to  be  "bought" 
with  increased  system  complexity. 

5.2.  COMBINATION  SYSTEM 

Experience  gained  in  studying  the  ternary  mode  pulse  torqued 
pendulous  integrating  accelerometer  leads  to  the  consideration  of  a 
combination  system.     Ternary  in  mode  of  error  detection,   that  is,   the 
system  senses  clockwise  error,   dead:  zone,  and  counter-clockwise 
error;    and  binary  in  mode  of  torquing,   that  is,   torque  clockwise  or 
counter-clockwise.     The  ternary  mode  in  error  detection  would  be 
used  to  effect  lead  in  the  torque  switching.     The  sense  of  torquing  would 
be  changed  on  entering  the  dead  zone.     Investigation  to  determine  if 
enough  lead  would  be  introduced  to  insure  lowest  order  limit  cycles  is 
recommended. 
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5.  3.     RECOMMENDATIONS  FOR  FURTHER  INVESTIGATION 

The  hardware  problems  specific  to  ternary  pulse  torquing  have 
not  been  investigated.     Investigation  to  determine  the  feasibility  of 
having  a  system  produce  a  reproducible  torque  pulse  in  two  senses  and 
also  not  produce  it  for  specific  periods  needs  to  be  carried  out.     An 
energy  sink  might  be  the  solution  except  in  the  case  when  the  ternary 
mode  is  used  to  measure  small  accelerations  at  little  expenditure  of 
energy.     The  environmental  effect  of  a  varying  rate  of  energy  expen- 
diture by  the  torquer  should  be  studied.     The  exact  performance  of 
error  detectors  at  high  sampling  rates  both  in  the  dead  zone  region  and 
away  from  the  dead  zone  is  needed  to  insure  the  feasibility  of  the  com- 
pensation schemes. 

Investigation  in  the  use  of  ternary  mode  pulse  torquing  gyros, 
gimbals,   and  other  inertial  instruments  is  recommended. 
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APPENDIX  A 
NOTATION 


AA  Dead  zone  of  error  detector 

A  Angular  Position  of  the  pendulum  with  respect  to  the  null 

position 

a  Component  of  acceleration  along  the  input  axis  of  the 

accelerometer 

aRA  Component  of  acceleration  along  the  Reference  Axis  (or  null 

axis)  of  the  accelerometer 

C  Coefficient  of  viscous  friction 

J  Moment  of  inertia  of  the  pendulum 

t  Time 

t  Time  between  samples 

M,  Torque  Generator  torque 

M  Magnitude  of  Torque  Generator  torque 

P  Pendulosity  of  the  Pendulum 

v  Velocity  increase 

v.  Indicated  velocity 

e  Total  error  in  velocity  indication 

e.  Digital  error 

e  Null  error 
n 

e  Cross-coupling  error 

6  Torquing  mode,   +1,   0,  or  -1. 

$  Non-dimensional  null  angle 

t  Non-dimensional  time 

<j)  d$/dr 

</>.  Non-dimensional  dead  zone 

r  Non-dimensional  acceleration  along  input  axis 

rR,  Non-dimensional  acceleration  along  reference  axis 

V  Non-dimensional  velocity 

E  Non-dimensional  velocity  error 

E  Non-dimensional  null  error 
n 

E  Non-dimensional  cross-coupling  error 
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E  ,  Non-dimensional  digital  error 

T  Non-dimensional  sampling  period 

<j>~  Freedom  in  position  of  a  simple  limit  cycle  about  a  switching  line 

s  Equal  by  definition 
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APPENDIX  B 

DERIVATION  OF  ERROR  DETECTOR  DEAD  ZONE  REQUIRED  TO 
PREVENT  A  LIMIT  CYCLE  WHEN  INPUT  ACCELERATION  IS  ZERO 

Consider  a  limit  cycle  of  K  time  intervals  in  the  positive  torque 
mode  followed  by  L  time  intervals  in  the  zero  torque  mode,  K  time  in- 
tervals in  the  negative  torque  mode  and  L  time  intervals  in  the  zero 
torque  mode.     The  limit  cycle  must  be  symmetric  in  that  the  sum  of  the 
torque  pulses  over  a  limit  cycle  must  be  zero.     Such  a  limit  cycle  is 
shown  in  Figure  B-l.     Two  requirements  must  be  met  in  order  that  the 
limit  cycle  can  exist: 

4<2K+L~1)  -*(K-1)=*A  (B_1) 

*(K+L-1)  "  *(2K+2L-1)  =  ^A  (B"2) 

Now  suppose  that  switching  line  1  is  held  fixed  and  switching  line  2 
is  moved  to  the  left  until  Condition  (B-l)  cannot  be  met.     Switching  will 
occur  prematurely  and  K  will  be  reduced  by  one.     Since  the  distance 
traveled  in  the  torque  mode  is  reduced,  the  distance  traveled  in  the  no 
torque  mode  must  increase  so  that  Condition  (B-l)  can  be  met  for  the 
new  value  of  K.     Thus  L  must  increase. 

Switching  line  2  can  be  moved  to  the  left  until  K  is  reduced  to  one 
and  L  is  increased  to  infinity.     As  L  approaches  infinity  the  velocity  at 
N  =  K  +  L  must  approach  zero  since  the  pendulum  has  been  damped 
wiftiout  being  forced  for  an  infinite  amount  of  time.     The  width  of  the 
limit  cycle  under  these  conditions  can  be  determined  from  Equations 
(2-7)  and  (2-8): 


"  *o  =  [J*  +  '"^  '  l     +    l  '  "J 


♦(CD)"*-    =       T-+e  6  I  (B- 


or 


♦(a,)  "  *o  ■  Ts  <B-4) 

Now,   if  switching  line  2  is  moved  to  the  left  such  that  </>.    is  greater 
than  T   ,   then  K  must  be  reduced  to  zero  and  the  limit  cycle  cannot  exist. 

-106- 


Positive 
Switch  Zone 


Negative 
Switch  Zone 


n  =  K  +  L 


n=  2K  +  L-1 


Figure  B-l.     Limit  Cycle  With  Zero  Acceleration. 
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Thus  the  condition  necessary  to  prevent  a  limit  cycle  is: 

*A   >  Ts  <B"5> 
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APPENDIX  C 

DERIVATION  OF  THE  FREEDOM  IN  POSITION  OF  A  SIMPLE 

LIMIT  CYCLE  ABOUT  THE  SWITCHING  LINE,  SHOWING   THE 

CONDITION  NECESSARY  FOR  A  SIMPLE  LIMIT  CYCLE  TO  EXIST 

A  simple  limit  cycle  is  defined  as  a  limit  cycle  which  is  com- 
plete after  two  switches.     For  the  purposes  of  this  derivation  it  is 
assumed  that  the  zero  switch  zone  is  of  sufficient  width  to  insure  that 
the  pendulum  in  steady  state  remains  in  the  positive  and  zero  switch 
zones  for  positive  acceleration  and  in  the  negative  and  zero  switch 
zones  for  negative  acceleration.     A  typical  limit  cycle  consisting  of 
K  sample  intervals  in  a  torque  mode  followed  by  L  sample  intervals 
in  the  zero  mode  is  shown  in  Figure  C-l. 

The  freedom  in  position  of  the  limit  cycle  is  defined  as  follows: 

*f  =    *(K  +  L-1)  "*<K-1)  (C'L) 

It  is  evident  from  Figure  C-l  that  the  complete  limit  cycle  may 
be  shifted  to  the  right  or  left  over  a  total  range  of  </>»  without  disturbing 
the  switching  sequence.     It  is  also  evident  that  if  </>f  is  negative,   the 
switching  sequence  must  be  altered  and  the  limit  cycle  described 
cannot  exist. 

The  velocity  summing  capability  of  the  system  depends  upon  the 
relation 

K 


(C.2) 


K  +  L 


Equation  (C.  2)  is  without  error  since  the  change  in  null  error  EN 
over  a  limit  cycle  is  zero.     Since  K  and  L  are  integers,    the  input  ac- 
celeration must  be  the  ratio  of  two  integers  if  the  limit  cycle  is  to  exist. 

The  equation  of  motion  of  the  pendulum  is: 

</>+</)=  (6  -  r)  (C.3) 

where  6  is  the  torquing  mode.     Equation  (C.  3)  is  linear  between  torque 
switches  and  may  be  integrated  as  such  with  results  being  valid  be- 
tween switches. 
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Positive 
Switch  Zone 


Zero 
Switch  Zone 


n=K 


n-K-1 


Negative 
Switch  Zone 


/*  n=K+L-l 


Figure  C-L     Simple  Limit  Cycle. 


-110- 


•  f-^^A  (     -{T-ri)) 


(C.4) 


</>T  =  0T    +</)T|l  -  e  x^+(6  -  r)(r  -  r.+e  -  lJ  (C.  5) 

A  solution  for  </>,„   T  *  in  terms  of  </>     is  obtained  using  Equation  (C.  4), 

'  .        -(K+L)T  -LT  -(K+L)T  /        -(K+L)T   \ 

(K+L)=*oe  +e  "e  -r(l-e  J      (C.6) 

where  T     is  the  non-dimensional  sampling  interval. 

Since  a  limit  cycle  is  complete  after  K  plus  L  sample  intervals, 

*o=*(K+L)  (C'7) 

Substituting  (C.  7)  into  (C.  6)  and  solving  for  </>    : 

-LT 

s       1 

4>     =  — ~-L~  +  (1  -  r)  (C.  8) 

°  -(K+DT 

1  -  e  s 

Substituting  (C.  2)  into  (C.  8) 

-LT 

p         s  -  1  T 

</>     = —  +  — ^-  (C.9) 

-(K+L)T  K  +  L 

-  e 

By  using  (C.  4),   (C.  5),   and  (C.  9)  an  appropriate  number  of  times 
-LTs       -(K+L-l)Ts  -(K-l)Ts  ^ 

*/K  n  =  4    + — ^-^ +  — ! (C 10) 

(K   l)         °  -(K+L)T  K  +  L 

1  -  e 

and 

-LTg       -(L-l)Tg  Tg  KT 

<W+T     n=4>n+- ~ ~l+e        +  (1  -  eT)  + —       (Cll) 

(K+L-l)         o  -(K+L)T  K  +  L 

1-  e  s 
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Substituting  (C.10)  and  (C.ll)  into  (C.  1) 

-(K+L-DT        -(K-l)T         -(L-l)T 

0f  -  1® 1-® --^ +  1  +  T    (. L  -  -*k_\        (C.  12) 

-(K+L)T  SV        K  +  L/ 

1  -  e  s  X  7 

Figures  C-2  through  C-6  show  </)f  as  a  function  of  the  non-dimensional 
sampling  period  for  various  values  of  K  and  L. 
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Figure  C-2.     Freedom  in  Position  of  Simple  Limit  Cycles  About  a 
Switching  Line;  1:N  Limit  Cycles. 
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Figure  C-3.     Freedom  in  Position  of  Simple  Limit  Cycles  About  a 
Switching  Line;    2:N  Limit  Cycles. 
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0.8  -' 


Figure  C-4.     Freedom  in  Position  of  Simple  Limit  Cycles  About  a 
Switching  Line;    3:N  Limit  Cycles. 
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Figure  C-5.     Freedom  in  Position  of  Simple  Limit  Cycles  About  a 
Switching  Line;    4:N  Limit  Cycles. 
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Figure  C-6.    Freedom  in  Position  of  Simple  Limit  Cycles  About  a 
Switching  Line;  5:N  Limit  Cycles. 
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APPENDIX  D 
SWITCHING  LOGIC   DESIGN 

D.  1.     INTRODUCTION  AND  DEFINITIONS 

The  object  of  this  appendix  is  to  present  the  relative  complexity 
of  the  different  switching  logics  required. 

The  logic  design  will  be  at  the  logical  element  level.     A  logical 
element  is  a  computer  element  which  performs  a  nontrivial  logic 
function;    that  is,   some  input  or  combination  of  inputs  produces 
something  new  as  an  output. 

In  this  presentation  of  switching  logic  design  elementary  know- 
ledge of  switching  design  is  presumed.     Chapter  17  of  Reference  6  is 
an  elementary  presentation  of  logical  design  and  was  used  in  de- 
veloping the  switching  logic  presented  here.     Switching  functions  will 
be  in  Boolean  algebra  form. 

Figure  D-l  defines  symbols  and  characteristics  of  the  logical 
elements  used.     The  delay  repeat  is  not  a  standard  element.     Figure 
D-2  is  an  elemental  layout  of  the  delay  repeat. 

Figure  D-3  shows  the  definition  of  the  error  level  signals.     The 
error  detector'  s  output  will  be  zero  in  the  dead  zone;   but  D,   the  dead 
zone  position  signal,   can  be  computed,   since  it  is  the  absence  of  the 
other  error  level  signals. 

Table  D-l  defines  the  letter  symbols  applicable  to  this  appendix 
only  used  to  define  quantities  used  in  the  logic  design. 

D.  2.     SWITCHING  LOGIC  FOR  ALTERNATING  PULSE  ZONE  COMPENSATION 

The  logic  for  switching  in  the  alternating  pulse  zone  compensation 
can  be  represented  by  Boolean  algebra  as  follows: 

Always  torque  when  in  the  torque  override  zones. 

P  =  F  (D.l) 

N  =  J  (D.  2) 

When  in  the  lead  compensation  zone,   torque  one  T   ,  do  not  torque 
next  T     if  still  in  the  lead  compensation  zone,  and  then  repeat  the  cycle 
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Figure  D-l.     Symbols  and  Characteristics  of  Logical  Elements, 
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Figure  D-2.     Elemental  Design  of  the 
Delay  Repeat. 
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TABLE  D-I 
DEFINITION  OF  LETTER  SYMBOLS  FOR  SIGNALS 
USED  IN  SWITCH  LOGIC    DERIVATIONS  ONLY 

A  -     Signal  A  at  present  time  instant 

A"         -     Not  A 

T+l 

A  -     Value  of  A  at  first  preceding  interval 

T+x 
A  -     Value  A  has  maintained  constant  for  x  preceding  periods. 

x  may  take  any  integer  value. 
P  -     Positive  torque 

N  -     Negative  torque 

D~        -     Error  in  dead  zone 

F  -     Error  in  positive  continuous  pulse  override  zone 

J  -    Error  in  negative  continuous  pulse  override  zone 

G  -     Error  in  positive  lead  compensation  zone 

H  -     Error  in  negative  lead  compensation  zone 

K  -     Positive  torque  pulse  as  a  result  of  G 

L  -     Negative  torque  pulse  as  a  result  of  H 

Q  -     Intermediate  signal  in  computing  K 

R  -     Intermediate  signal  in  computing  L 

(+or  -)  _      -     Mean  error  bias  correction  magnitude  when  switching  about 

</>  equal  (-  or  +)  </>„ 
(+or-)  V      -     Mean  error  bias  correction  magnitude  less  U/2  when  switching 

about  </>  equal  (-  or  +)  </>T 

W         -     Negative  error  sense 

X  -     Permit  positive  torquing  for  more  than  one  T 

X  -     Positive  torque  for  not  more  than  one  T 

^  s 

Y  -     Permit  negative  torquing  for  more  than  one  T 

Y  -     Negative  torque  for  not  more  than  one  T 

B  -     Error  has  been  in  the  dead  zone  this  sampling  instant 

and  the  one  previous 
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as  long  as  in  the  lead  compensation  zone. 


K  =  G  •   KT+1  (D.  3) 


Combining 


L  =  H  •    LT+1  (D.4) 


P  =  F  +  K  (D.4) 


P  =  F  +G  •   KT+1  (D.  5) 


N  =  J  +  L  (D.6) 


N  =  J  +  H  ♦   LT+1  (D.7) 

The  switching  logic  described  by  equations  (D.  5)  and  (D.7)  is 
shown  in  Figure  D-4. 

D.  3.     SWITCHING  LOGIC  DESIGN  FOR  INVERSE  ORDER  ALTER- 
NATING PULSE  ZONE  COMPENSATION 

The  switching  logic  for  the  inverse  order  alternating  pulse  zone 
compensation  is  like  the  switching  logic  for  the  alternating  pulse  zone 
compensation  except  that  when  in  the  lead  compensation  zone  the 
opposite  action  of  that  taken  the  preceding  period,  no  matter  which 
zone  determined  the  action,   is  desired.     This  is  done  by  changing 
equation  (D.  3)  to 


K  =  G  •    PT+1  (D.  8) 


and  equation  (D.  4)  to 


Then 


T  +1 
L  =  H  •   N       .  .  (D.  9) 


P  =  F  +  G  •    PT+1  (D.  10) 


N  =  J  +H  •   NT+1  (D.ll) 

The  switching  logic  described  by  equations  (D.  10)  and  (D.  11) 
is  shown  in  Figure  D-5. 
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Figure  D-4.     Switching  Logic  For  Alternating  Pulse  Zone  Compensation. 
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Figure  D-5.     Switching  Logic  for  Inverse  Order  Alternating 
Pulse  Zone  Compensation. 
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D.  4.     LOGICALLY  FORCED  LOWEST  ORDER  COMPENSATION 
SWITCHING  LOGIC 

In  the  logically  forced  lowest  order  compensation,   for  values 
of  r  less  than  0.  5  the  system  is  allowed  to  torque  for  only  one  sam- 
ple period,   T   ,   at  a  time,  and  for  values  of  r  greater  than  0.  5  the 

fa 

system  is  allowed  to  not  torque  only  one  sample  period,   T   ,   at  a 

fa 

time.     The  logic  will  be  changed  from  multiple  T     positive  torquing, 

s 

X,    to  single  T     positive  torquing,   X;   or  from  multiple  T     negative 

fa  fa 

torquing,   Y,    to  single  T     negative  torquing,   Y;     when  the  error 
signal  remains  in  the  dead  zone  for  two  consecutive  sample  periods. 
This  is  represented  by  Booloean  algebra  by 

B    =   D  •    DT+1  ,  (D.  12) 


reset  X 


and  reset  Y 


X   =  BT+X  ,  (D.13) 


Y    =  BT+X  (D.  14) 

The  logic  will  be  changed  from  single  T     torquing,  X  or  Y, 

fa 

to  multiple  T     torquing,  X  or  Y  respectively,   when  the  error  signal 
passes  the  -h/>,    or  -</>_    line  respectively  in  an  increasing  error 
direction. 


Set  X 
g 


Set  Y„ 
g 


X  =  FT+X.  (D.  15) 


Y  =  JT+X  (D.  16) 


The  torque  switching  can  be  summarized  as  follows: 

—     — T+l 
P  =  <F  +  G  •   X  +  G  •   X  •    P    r  (D.  17) 

N  =   J  +  H  •    Y+H-T-   NT+1  (D.  18) 
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The  logic  described  by  equations  (D.  12)  through  (D.18)  for  the 
logically  forced  lowest  mode  compensation  is  shown  in  Figure  D-6. 

D.  5.     ERROR  BIAS  CORRECTION  SWITCHING  LOGIC 

Proceeding  from  positive  error  to  negative  error,   the  change 
in  error  bias  correction  is  +V  when  changing  from  +<L    to  +<j>-n, 
+U  when  changing  from  +</>„  to  -</»n,   and  +V  when  changing  from 
-<£      to  -<£    .     The  sign  of  the  change  in  error  bias  correction  is 
negative  proceeding  from  negative  to  positive.     The  magnitude  of 
U  and  V  may  be  adjusted  to  make  the  best  correction  for  the  mean 
error  due  to  switching  about  a  line. 

When  it  is  necessary  to  correct  for  the  error  only  when  switching 
about  the  <£,    lines,   the  Boolean  algebra  equations  for  the  switching  are 

fV+^1    =  F-    DT+X+D-    LT+X  (D.19) 


/  V  +  ^  1   =  D  •    FT+X  +  L  •  DT+X  (D.  20) 


Only  one  quantity  may  be  stored  at  a  time.     Figure  D-7  shows  the 
required  switching  logic. 

If  it  is  necessary  to  correct  only  for  the  error  bias  in 
switching  about  the  </>_  lines  as  in  the  uncompensated  case  or  the 
logically  forced  lowest  mode  case  again  a  special  starting  condition 
need  not  be  used  unless  the  pendulum  was  restrained  until  counting 
started.     This  is  because  the  pendulum  will  be  switching  about  the 
initial  switch  line  when  counting  starts.     +U  must  be  switched  when 
entering  the  F  or  G  zones  when  -U  was  last  switched.     -U  must  be 
switched  when  entering  the  H  or  L  zones  when  +U  was  last  switched. 

+U  =  (F  +  G)  •  (-U)T+X  (D.  21) 


-U  =  (H  +  L)  •    (+U)T+X  (D.  22) 
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Figure  D-6.     Torque  Switching  Logic  for  Logically  Forced 
Lowest  Order. 
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Figure  D-7.     cj)^  Switching  Line  Error  Bias  Correction 
Switching  Logic. 
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-U  will  reset  (+U)T+X  and  +U  will  reset  (-U)T+X. 

Figure  D-8  shows  the  two  state  </>~  switching  line  error  bias 
correction  switching  logic. 

A  four  state  error  bias  correction  to  correct  for  the  error 
bias  due  to  switching  about  each  of  four  possible  switching  lines  can 
be  constructed     using  a  combination  of  the  two  preceding  error  bias 
corrections. 
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Figure  D-8.     </»D  Switching  Line  Error  Bias  Correction 
Switching  Logic. 
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